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En el presente trabajo se desarrolla el diseño de un 
controlador robusto de orden fraccional —el cual se denota 
como FPI, debido a que constituye una generalización del 
controlador PI convencional— para el control efectivo de 
la distribución de agua en canales principales de riego 
caracterizados por presentar una amplia variación en 
sus parámetros dinámicos. El interés en esta clase de 
controladores se justifica por hecho de que la dinámica 
de los canales principales de riego puede cambiar de 
forma drástica con las variaciones de sus regímenes de 
descargas. El controlador FPI desarrollado se programó 
en un PLC SIMATIC S-7 300 y se implementó en tiempo 
real en el primer tramo del canal principal Imperial de 
Aragón, perteneciente a la Confederación Hidrográfica del 
Ebro, España. Se realizaron experimentos comparativos 
de funcionamiento en tiempo real del controlador FPI 
diseñado con otro controlador PI equivalente, en el sentido 
de exhibir ambos el mismo comportamiento dinámico 
en lazo cerrado para las especificaciones temporales de 
la planta nominal. Los resultados obtenidos de estos 
experimentos demuestran que cuando los parámetros 
dinámicos del tramo del canal presentan una amplia 
variación, el controlador FPI exhibe una mayor robustez, 
lo cual posibilita lograr una mayor velocidad de respuesta, 
satisfacer de forma operativa las demandas de los usuarios 
y minimizar las pérdidas por concepto de operación de los 
recursos hidráulicos disponibles. 


Palabras clave: controlador robusto de orden fraccional, 
implementación práctica en tiempo real, canal principal de 
riego, parámetros dinámicos variables en el tiempo, uso 
eficiente de recursos hidráulicos. 


Abstract 


RIVAS-PÉREZ, R., FELIU-BATLLE, V., CASTILLO-GARCÍA, 
E, SÁNCHEZ-RODRÍGUEZ, L. €s LINARES-SÁEZ, A. Robust 
fractional order controller implemented in the first pool of the 
Imperial de Aragón main canal. Water Technology and Sciences 
(in Spanish). Vol. V, No. 1, January-February, 2014, pp. 23-42. 


This paper presents the design of a robust fractional order 
controller —which we denote as EPI controller since it can be 
regarded as a generalization of the standard Pl controller— for 
effective control of water distribution in irrigation main canal 
pools characterized by time-varying dynamical parameters. 
Interest in such robust fractional order controller is justified by 
the fact that dynamical parameters of irrigation main canal pools 
may change drastically according to their operating discharge 
regimes. The FPI controller designed was installed in a PLC 
SIMATIC S-7 300, and was real-time implemented in the first 
pool of the Imperial de Aragon main canal, belonging to the Ebro 
Hydrographical Confederation (Spain). Real-time experiments 
were carried out with equivalent —in the sense of exhibiting the 
same closed loop dynamics for the nominal plant specifications— 
FPI and Pl controllers. They showed a better performance and 
robustness of the FPI controller when large canal parameters 
variations were produced, that which facilitates to achieve a 
bigger speed of response, to satisfy in an operative way the 
demands of the users, as well as to minimize the losses for concept 
of operation of the available hydraulic resources. 

Keywords: robust fractional order controller, real-time 
implementation, irrigation main canal, large time-varying 
dynamical parameters, efficient management of hydraulic 
resources. 
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Introducción 


La experiencia mundial muestra que una vía 
importante para aumentar la operatividad 
y eficiencia en la distribución de agua en los 
canales principales de riego consiste en la 
aplicación de sistemas efectivos de control 
automático y de ayuda a la toma de decisión, 
los cuales posibilitan la obtención de una 
estricta correspondencia entre las demandas 
de los diferentes usuarios y las entregas de 
este recurso, así como reducir las pérdidas de 
agua por concepto de operación hasta en un 
90% (Clemmens, 2006; Feliu-Batlle ef al., 2007; 
Kovalenko, 1983; Litrico y Fromion, 2009; 
Rivas-Pérez et al., 2003). 

El objetivo fundamental del desarrollo e 
implementación práctica de sistemas efecti- 
vos de control automático de la distribución 
de agua en los canales principales de 
riego consiste en satisfacer, a pesar de las 
incertidumbres, las demandas de agua de 
los diferentes usuarios, garantizando que los 
volúmenes de agua extraídos desde la fuente 
de abasto se correspondan con las necesidades 
reales, así como que se minimicen las pérdidas 
por concepto de operación en todo el canal 
(Buyalski et al., 1991; Kovalenko, 1983; Rivas- 
Pérez, 1990). 

Las investigaciones desarrolladas en dife- 
rentes canales principales de riego muestran 
que sus parámetros dinámicos experimentan 
grandes variaciones cuando sus regímenes 
de descarga varían en el rango de operación 
[O Q...] y/o otros parámetros hidráulicos 
cambian (Corriga et al., 1989; Litrico y Fromion, 
2009; Montazar et al., 2005; Rivas-Pérez et al., 
2007). Esta clase de canales se conoce como 
canales con parámetros dinámicos variables 
en el tiempo (Litrico et al., 2006; Rivas-Pérez 
et al., 2008b). Por consiguiente, cualquier 
controlador que se diseñe debe comportarse de 
forma robusta frente a esta clase de variaciones 
(Deltour y Sanfilippo, 1998; Feliu-Batlle et al., 
2009a; Litrico et al., 2006). 

Los controladores PID (proporcional-inte- 
gral-derivativo) constituyen los controladores 


más utilizados y expandidos en los sistemas 
de control automático de la distribución 
de agua en los canales principales de 
riego (Buyalski et al., 1991; Clemmens y 
Schuurmans, 2004; Litrico y Fromion, 2009; 
Malaterre et al., 1998; Schuurmans et al., 1999). 
Sin embargo, diferentes estudios muestran 
que cuando los canales principales de riego 
presentan parámetros dinámicos variables 
en el tiempo, la aplicación de controladores 
PID 
originándose un aumento significativo en el 


convencionales no resulta efectiva, 
tiempo de establecimiento y en el pico máximo 
de la respuesta temporal del sistema de control, 
lo cual implica que disminuya la operatividad 
de la distribución de agua en los canales y, por 
consiguiente, que aumenten las pérdidas de los 
recursos hidráulicos disponibles por concepto 
de operación (Clemmens y Schuurmans, 2004; 
Litrico y Fromion, 2009; Feliu-Batlle et al., 
2009b; Rivas-Pérez et al., 2002; Wahlin, 2004). 
Es por ello que el desarrollo de controladores 
robustos de la distribución de agua en canales 
principales de riego constituye un área de 
investigación de elevado interés científico- 
técnico y práctico (Litrico y Fromion, 2009; 
Wahlin y Clemmens, 2006). El diseño de esta 
clase de controladores constituye un reto y un 
problema vital para mejorar la operatividad 
y eficiencia en la distribución de agua en los 
canales principales de riego. 

En los últimos años, como resultado de 
una mejor comprensión del cálculo fraccional 
y de la salida al mercado de un nuevo 
circuito electrónico conocido como fractal, 
los operadores de orden fraccional han 
sido aplicados con resultados satisfactorios 
en el modelado y control de procesos con 
comportamiento dinámico complejo, entre los 
que se encuentran los procesos con parámetros 
distribuidos (Chen et al., 2004; Machado, 1997; 
Petras, 2002; Podlubny, 1999a; Vinagre et al., 
2000). El cálculo fraccional representa al campo 
de la matemática que involucra a derivadas 
e integrales de orden no entero (arbitrario), y 
constituye una generalización de los conceptos 
estándar de diferenciación e 
(Podlubny, 1999a). 


integración 
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La aplicación de controladores de orden 
fraccional (FOC) posibilita ampliar las acciones 
de control que se pueden desarrollar sobre una 
variada clase de plantas y /o procesos, así como 
resolver el problema de la automatización 
efectiva de procesos, cuyo comportamiento 
dinámico se describe mediante modelos de 
orden no entero, sin utilizar aproximaciones 
(Podlubny, 1999a). Una de las ventajas de 
esta clase de controladores consiste en que 
posibilitan el diseño de sistemas de control 
robusto de plantas, cuyos parámetros dinámicos 
presentan un amplio rango de variación o de 
plantas sometidas a grandes perturbaciones 
(variaciones en la carga), así como de plantas 
con grandes retardos de tiempo, logrando 
mantener la estabilidad del sistema de control 
y presentando una simple implementación 
práctica (Feliu-Batlle et al., 2009b; Monje et al., 
2004; Vinagre et al., 2000). 

El comportamiento cualitativo de los 
controladores PID puede ser mejorado de 
forma considerable mediante su generalización 
como controlador de orden fraccional PI*D”, 
que involucra un integrador de orden a y un 
diferenciador de orden A (Podlubny, 1999b). En 
consecuencia, los controladores PI*D* han sido 
propuestos y han recibido una considerable 
atención (Barbosa et al., 2004; Castillo et al., 
2010; Chen et al., 2006; Djouambi et al., 2005; 
Leu et al., 2002; Monje et al., 2004; Petras et al., 
2003; Podlubny et al., 2002; Rivas-Pérez et al., 
2012; Wang et al., 2005). 

De manera reciente han sido reportados 
trabajos relacionados con la aplicación de 
controladores de orden fraccional en el control 
de la distribución de agua en canales principa- 
les de riego (Castillo-García et al., 2013; Feliu- 
Batlle et al., 2011, 2009a, 2009b, 2007; Monje 
et al., 2004; Rivas-Pérez et al., 2008c; Sánchez 
et al., 2007). Sin embargo, los resultados que 
se ofrecen en estos trabajos presentan en 
general un carácter teórico y no se encuentran 
implementados en canales principales de riego 
en explotación. 

Diferentes autores han reportado que a 
pesar de que el problema del diseño de sis- 


temas efectivos de control de la distribución 
de agua en canales principales de riego ha sido 
objeto de numerosas publicaciones científicas, 
sólo una pequeña parte de estos controladores 
ha sido implementada de forma exitosa en la 
actividad práctica (Clemmens y Schuurmans, 
2004; Litrico y Fromion, 2009; Schuurmans et 
al., 1999; Wahlin y Clemmens, 2006). 

El controlador PI (proporcional integral) es 
el más utilizado en el control de la distribución 
de agua en los canales principales de riego 
debido a que es más fácil de ajustar de 
forma correcta que un PID y a que la parte 
derivativa del PID es muy sensible al ruido 
de los sensores, originando problemas de 
estabilidad del sistema de control (Baume et 
al., 1999; Burt et al., 1998; Rivas-Pérez et al., 
1998). El método habitualmente utilizado para 
ajustar los parámetros de estos controladores 
es el de prueba y error (Deltour y Sanfilippo, 
1998; Seatsu, 1999; Whalin y Clemmens, 
2006), o mediante optimización (Clemmens 
y Schuurmans, 2004). Estos métodos, por lo 
general, se basan en modelos matemáticos 
de la planta nominal, y no consideran que 
los parámetros dinámicos de los canales 
principales de riego pueden cambiar en 
un amplio rango, como resultado de las 
variaciones de sus regímenes de descarga y/o 
de otros parámetros hidráulicos, por lo que los 
controladores PI se desajustan con facilidad o 
incluso pueden llegar a inestabilizar el sistema 
de control (Litrico et al., 2006; Rivas-Pérez 
et al., 2007). Es por ello que en este trabajo se 
exploran las posibles ventajas de un regulador 
de orden fraccional FPI versus un regulador 
convencional PI desde el punto de vista de su 
robustez frente a variaciones de los parámetros 
dinámicos del canal y de las perturbaciones. 

Este trabajo tiene como objetivos: a) el 
desarrollo de un controlador robusto de orden 
fraccional FPI para el control efectivo de la 
distribución de agua en un tramo de un canal 
principal de riego (caracterizado por presen- 
tar grandes variaciones en sus parámetros 
dinámicos), que posibilite disminuir el tiempo 
de establecimiento y el pico máximo de la 
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respuesta temporal del sistema de control, 
lo cual implica aumentar la operatividad de 
la distribución de agua y, por consiguiente, 
reducir las pérdidas por concepto de operación 
de los recursos hidráulicos disponibles; b) 
la implementación y validación práctica del 
controlador FPI diseñado en el primer tramo 
del canal principal Imperial de Aragón. 

En este trabajo se considera que es posible 
desarrollar sistemas efectivos de control des- 
centralizado de la distribución de agua en 
canales principales de riego mediante el 
diseño de controladores robustos de orden 
fraccional independientes para cada tramo 
del canal y adicionar desacopladores, con el 
objeto de minimizar la interacción entre los 
tramos adyacentes (Litrico y Fromion, 2009; 
Schuurmans et al., 1999; Wahlin y Clemmens, 
2006). Es por ello que este trabajo se centra 
en el diseño de un controlador FPI para un 
solo tramo de un canal principal de riego. El 
problema del control efectivo de la distribución 
de agua en tramos adyacentes o en canales 
completos (íntegros) mediante controladores 
de orden fraccional se analizará en trabajos 
futuros de investigación. 

La fundamental contribución de este tra- 
bajo consiste en la implementación y valida- 
ción práctica (por primera vez) de un con- 
trolador robusto de orden fraccional (FPI) de 
la distribución de agua en un tramo de un 
canal principal de riego en explotación (primer 
tramo del canal principal Imperial de Aragón) 
con resultados muy satisfactorios. Esta con- 
tribución constituye una solución práctica y 
relativamente simple al complejo problema 
del diseño de controladores robustos de la 
distribución de agua en canales principales de 
riego con comportamiento dinámico impreciso 
(con grandes incertidumbres), y verifica que 
los controladores de orden fraccional (FOC) 
posibilitan obtener mejores resultados que los 
controladores convencionales (PI) en esta clase 
de aplicación. 

Los resultados obtenidos formaron parte 
de un proyecto de investigación entre la 
Universidad de Castilla-La Mancha y la 


Confederación Hidrográfica del Ebro (España) 
en modelado y control de la distribución de 
agua en el canal principal Imperial de Aragón. 

El trabajo se encuentra organizado de 
la siguiente forma: se obtiene un modelo 
matemático orientado al control del compor- 
tamiento dinámico del primer tramo del canal 
principal Imperial de Aragón; se desarrolla 
el diseño en el dominio de la frecuencia de 
un controlador robusto de orden fraccional 
(FPI) y de un controlador PI bajo las mismas 
especificaciones de respuesta temporal para 
el tramo de canal objeto de estudio bajo 
régimen de operación nominal; se analizan 
los resultados experimentales obtenidos 
de la implementación práctica de ambos 
controladores (FPI y PI) en el primer tramo 
del canal principal Imperial de Aragón, y 
finalmente se presentan las conclusiones del 
trabajo. 


Modelo matemático orientado al control 
de un tramo de un canal principal de 
riego 


La necesidad de disponer de modelos ma- 
temáticos que describan de forma adecuada 
el comportamiento dinámico más relevante 
de los canales principales de riego constituye 
uno de los inconvenientes que presenta el 
diseño e implementación práctica de sistemas 
efectivos de control de la distribución de agua 
en los canales principales de riego (Litrico y 
Fromion, 2009; Kovalenko, 1983; Rivas-Pérez 
et al., 2008a). Esta clase de modelos ofrece una 
valiosa información sobre las propiedades de 
los sistemas dinámicos, tales como respuesta 
temporal, retardo de tiempo, constantes de 
tiempo, frecuencia de cruce, pico máximo, 
etcétera. 

El modelado matemático del compor- 
tamiento dinámico de los canales principales 
de riego se realiza usualmente mediante las 
ecuaciones de Saint-Venant, las cuales se 
basan en las leyes clásicas de conservación de 
masa y de conservación del momento (Chow, 
1988; Kovalenko, 1983; Rivas-Pérez et al., 
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2003), y representan ecuaciones no lineales 
en derivadas parciales del tipo hiperbólicas 
con restricciones complejas, por lo que su 
utilización de forma directa en el diseño de 
sistemas de control presenta serias dificultades 
(Litrico y Fromion, 2009). Esto fundamenta la 
necesidad de obtener modelos matemáticos 
simples de los canales principales de riego que 
posibiliten su aplicación directa en el diseño de 
sistemas efectivos de control. 

son usualmente 
suficientes para capturar las propiedades diná- 


Los modelos lineales 


micas fundamentales de los canales principales 
de riego, requeridas en el diseño de sistemas de 
control (Litrico y Fromion, 2009; Rivas-Pérez et 
al., 2007). Esta clase de modelos puede obtenerse 
mediante la aplicación de las herramientas de 
identificación de sistemas (Rivas-Pérez et al., 
2011). La identificación de sistemas permite 
obtener modelos matemáticos simples que 
describen de forma adecuada la dinámica de 
plantas complejas para su utilización en el 
diseño de sistemas efectivos de control (Ljung, 
1999% 

Un canal principal típico de riego se 
encuentra conformado por diferentes tramos 


separados por compuertas, las cuales son 
utilizadas para regular la distribución de agua 
de un tramo a otro (ver figura 1, en la cual 
se muestra un diagrama del primer tramo 
de canal principal Imperial de Aragón). La 
magnitud de apertura de las compuertas se 
varía de forma adecuada para mantener un 
perfil determinado del agua a través de todo 
el tramo. En general, en un canal principal 
de riego controlado automáticamente (me- 
diante el método de regulación aguas abajo 
con sensor de nivel alejado) la variable 
controlada es el nivel de agua y(t) medida 
al final del tramo, la variable de control es la 
posición de la compuerta u(t) aguas arriba 
y las perturbaciones fundamentales son las 
descargas laterales q (t), donde ¡=1, 2, ..., n es el 
número del tramo del canal (Kovalenko, 1983; 
Malatertre et al., 1998). 

Los datos y resultados reportados en este 
trabajo se obtuvieron del tramo 1 (origen) del 
canal principal Imperial de Aragón, conocido 
como Bocal. Este canal forma parte de la 
Confederación Hidrográfica del Ebro, la cual 
constituye la confederación hidrográfica más 
importante de la península Ibérica (Pedregal 


Yr (t) 
Compuerta 
equivalente 


Yi (0) 


LS = sensor de nivel. 
GPS = sensor de posición de compuerta. 


Figura 1. Diagrama del tramo 1 (bocal) del canal principal Imperial de Aragón. 
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et al., 2009). Este canal toma sus aguas por 
derivación del río Ebro. 

El Bocal presenta una longitud de 8 kiló- 
metros, una profundidad variable entre 3.1 
y 3.6 metros, un ancho variable entre 15 y 
26.9 metros, y conduce un caudal máximo 
(descarga) de Q(t) = 30 m*/s en toda su 
longitud. En la figura 2 se ofrece una vista 
superior de un segmento del Bocal, en la cual 
es posible observar los grandes volúmenes de 
agua que se transportan a través de este tramo 
de canal. 

El nivel de agua en el Bocal se controla a 
través de 10 compuertas sumergibles mediante 
el método de regulación aguas abajo con 
sensor de nivel alejado. Estas compuertas se 
encuentran ubicadas en un recinto denominado 
como Casa de Compuertas. Las medidas 
disponibles en el Bocal son nivel de agua aguas 
arriba (río Ebro) y,(t), nivel de agua aguas 
abajo y,,(t), nivel de agua aguas abajo alejado 
y,(t), y posición (abertura) de las compuertas 
aguas arriba u,(t) (ver figura 1). En la figura 
1 se representan las 10 compuertas situadas 
en la Casa de Compuertas mediante una sola 
compuerta equivalente (Rivas-Pérez et al, 
2008a). 

Con el objeto de obtener el modelo 
matemático que describe el comportamiento 
dinámico de la distribución de agua en el Bocal 
(planta) se desarrollaron experimentos basados 


Figura 2. Vista superior de una porción del bocal. 


en la respuesta escalón. Para ello, se mantuvo 
en posición fija la compuerta aguas abajo; 
un total de cuatro compuertas aguas arriba 
recibieron un incremento simultáneo en su 
magnitud de apertura de 25 cm, lo que equivale 
a un incremento en la magnitud de apertura 
total de 100 cm, y se midió la variación de 
nivel de agua aguas abajo alejado. La respuesta 
temporal obtenida del comportamiento 
dinámico del Bocal se presenta en la figura 3. 

Esta respuesta muestra que el compor- 
tamiento dinámico de la distribución de agua 
en el Bocal puede representarse mediante la 
siguiente función de transferencia de segundo 
orden con retardo de tiempo: 


Ñ Ay (5) Ñ K —TS 
EA Au (s) (s+1(Bs+1) 1 


donde Ay,,(s) es la variación de nivel de agua 
aguas abajo alejado; Au,(s), la variación de 
la magnitud de apertura de la compuerta 
equivalente; K, la ganancia estática; T,, T,, 
las constantes de tiempo, y 1 es el retardo de 
tiempo. Se considera que T, es la constante de 
tiempo dominante (la de mayor magnitud, 
asociada con la dinámica del tramo del canal), 
mientras que T, es la constante de tiempo de 
menor magnitud, la cual representa la dinámica 
secundaria (motor + compuerta + sensor de 
nivel). Usualmente, T, es mucho menor que T, 
(Rivas-Pérez et al., 2007). 

Cuando la descarga Q,(t) a través de 
las compuertas aguas arriba del Bocal se 
corresponde con el régimen nominal de ope- 
ración Q,(t) = Q,,,, se obtienen los valores 
nominales de los parámetros del modelo 
matemático (1) (planta nominal). Estos pará- 
metros se representan como K,, Ty Tay Ty 
y sus valores son los siguientes K, = 0.055, 
T¡, = 880.79 s, T,, = 81.27 s, 1, = 360 s. Los 
resultados de validación del modelo lineal 
nominal obtenido se muestran en la figura 4. 
De esta figura se observa que existe una buena 
correspondencia entre la salida medida y las 
predicciones que realiza dicho modelo. 
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Centímetros 


Centímetros 


Bocal, respuesta temporal (variación de nivel aguas abajo alejado) 


¡OS 


0 1000 2 000 3 000 4 000 5 000 6 000 
Tiempo (s) 
Bocal, variación apertura compuerta 


0 1 000 2 000 3000 4 000 5 000 6000 
Tiempo (s) 


Figura 3. Respuesta temporal del bocal frente a una señal escalón en sus compuertas aguas arriba. 


Bocal, validación modelo 


Salida medida 
== == Salida del modelo 
Índice de ajuste: 93.79% 


Variación nivel de agua (cm) 


0 1000 2 000 3 000 4 000 5 000 6 000 
Tiempo (s) 


Figura 4. Resultados de validación del modelo matemático nominal de la distribución de agua en el bocal. 
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Sin embargo, cuando la descarga Q,(t) a 
través de las compuertas aguas arriba del Bocal 
no se corresponde con el régimen nominal de 
operación Q.(t) + Q_.. y presenta variaciones 
en el rango de operación [Q,., Q..] se 
originan cambios: en el volumen de agua 
que se acumula en el tramo del canal entre el 
volumen de descarga mínima y el de máxima 
descarga; en el nivel de agua aguas abajo en el 
rango [H,., H,...], así como en la propagación 
del flujo, la cual se determina en función de la 
celeridad de la ola (C) y de la velocidad media 
del flujo (V) (Chow, 1988; Kovalenko, 1983). 
Como resultado, los parámetros dinámicos 


mín 


del Bocal experimentan una amplia variación 
(incertidumbres paramétricas del modelo) en 
los siguientes rangos: 


0.01 <K t) <0.16 


600 <T, (+) <3000 


2) 


50 s T, (+) s 200 


) 
300 s 1(t) < 360 


Por esta razón, cualquier controlador que se 
diseñe para el Bocal debe garantizar a priori un 
nivel específico de comportamiento mínimo en 
todo el rango de variación de sus parámetros 
dinámicos (incertidumbres del modelo). Este 
es el problema del diseño de controladores con 
comportamiento robusto. 


Diseño del controlador robusto de orden 
fraccional (FPID) 


En esta sección se desarrolla el diseño en el 
dominio de la frecuencia de un controlador 
robusto deorden fraccional FPIdeladistribución 
de agua para el tramo del canal principal, 
objeto de estudio bajo comportamiento no- 
minal. Como se mencionó previamente, este 
controlador puede ser considerado como una 
generalización de los controladores PI. Por otra 
parte, se analizan las propiedades de robustez 
del controlador FPI en comparación con las de 
un controlador convencional PI equivalente. 


En el Bocal se encuentra implementado 
un controlador PI, el cual presenta como 
ventajas su facilidad de ajuste (sintonización) y 
obtención de un error en régimen permanente 
nulo frente a cambios tipo escalón en la 
referencia O perturbaciones también en forma 
de escalón. Como se señaló previamente, estos 
controladores presentan la dificultad de ser muy 
sensibles a las variaciones de los parámetros 
dinámicos del canal, desajustándose con 
facilidad oO incluso pudiendo llegar a 
inestabilizar el sistema de control (Litrico y 
Fromiom, 2009; Rivas-Pérez, 1990). 

El enfoque desarrollado en este trabajo 
consiste en diseñar un controlador que exhiba 
el mismo comportamiento que un PI frente 
al comportamiento dinámico nominal del 
Bocal; es decir, diseñado para las mismas 
especificaciones, pero que presente menos 
sensibilidad frente a cambios en los parámetros 
dinámicos de la planta y, por consiguiente, 
mayor robustez. Nuestra propuesta novedosa 
consiste en la utilización de un controlador de 
orden fraccional FPI, donde el grado adicional 
de libertad que posibilita aumentar la robustez 
y conseguir comportamientos frecuenciales 
prefijados viene dado por el orden fraccional 
a del controlador. Por ello, nuestro objetivo 
consiste en obtener un controlador que, 
comportándose de forma similar a un PI bajo 
régimen nominal, proporcione una mejor 
respuesta temporal del sistema de control 
cuando se produzcan variaciones en el régimen 
de descarga a través de las compuertas aguas 
arriba del Bocal en el rango de operación [O ;., 
Qu.) y/o variaciones de otros parámetros 
hidráulicos. 

El controlador FPI que se propone se 
representa mediante la siguiente función de 


transferencia: 
al 
Repj(s)=K, HL (3) 
donde: 
a no 
gs (jo) =0we 2 (4) 
0O<as<2 
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Observar que el controlador FPI (3) se 
convierte en un controlador PI convencional 
cuando el parámetro de orden fraccional a = 
1. Esta clase de controladores requiere el ajuste 
de tres parámetros K,, T, y a, es decir, un 
parámetro más que en el caso del controlador PI 
convencional. El parámetro de orden fraccional 
a se utiliza para cumplir con especificaciones 
adicionales de comportamiento robusto del 
sistema de control. 

En la figura 5 se muestra el diagrama 
de bloques del sistema de control de orden 
fraccional de la distribución de agua en el 
Bocal que se propone. En este diagrama se 
ha considerado una perturbación D(s), la cual 
representa la descarga lateral del tramo q,(t) 
(ver figura 1) y se modela mediante una señal 
escalón que pasa a través de un filtro de primer 
orden con constante de tiempo T, =70 s. 

El método que se propone de diseño 
del controlador robusto FPI se basa en las 
siguientes especificaciones en el dominio 
de la frecuencia: a) margen de fase deseado 
(4, ), el cual garantiza el amortiguamiento 
nominal deseado de la respuesta temporal 
(sobreimpulso (M),)), así como la robustez 
del sistema de control frente a variaciones 
del retardo de tiempo; b) frecuencia de cruce 
de ganancia deseada (w,), la cual garantiza la 


velocidad nominal deseada de la respuesta 
temporal del sistema de control (tiempo de 
establecimiento (t)), c) zero error de estado 
estacionario. La última especificación implica 
que el controlador debe incluir un término 
integral. Los parámetros del controlador (3) 
que cumplen con las especificaciones a)-c) se 
calculan mediante el siguiente procedimiento 
algebraico desarrollado en el dominio de la 
frecuencia. 

En correspondencia con el teorema del valor 
final (Ogata, 1993), la especificación c) implica 
que a > 0. Las especificaciones de diseño a) y 
b) dadas por un determinado margen de fase 
(4, ) y una frecuencia de cruce de ganancia 
(w.) pueden expresarse de forma compacta 
utilizando la siguiente notación compleja: 


Repr (jo.)G (jo) E e (-Am)i = ión (5) 


donde G(¡w.) representa al comportamiento 
dinámico (1) del Bocal. 

De (5), considerando las expresiones (3) y 
(4), se obtiene: 


TU 
¡Za : 
: 1+T,0%e 2 —elón 
Repi(j0.)=K, md == (6) 


Perturbación tipo 
señal escalón 


O) 


Compuerta + dinámica del canal 


(T,s+1XT,s+1) 


Figura 5. Diagrama de bloques del sistema de control de orden fraccional de la distribución de agua en el bocal. 
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Teniendo en cuenta que: 


(jo,)"=0w2 [cos[5a + sin(£a) (7) 


de la expresión (6) se obtienen los parámetros 
del controlador K, y T; 


donde: 
. _ R Joel (10) 
G(jo,) 
e (11) 
G(jo,) 
RO y SO representan las partes real 


e ¡imaginaria de un número complejo, 

respectivamente. 

Las especificaciones de diseño del 
controlador FPI son las siguientes: 

e Margen de ganancia: el parámetro 
que presenta la mayor influencia en la 
estabilidad del sistema de control en lazo 
cerrado es la ganancia estática K de la 
planta. Los resultados presentados en la 
sección anterior (2) muestran que este 
parámetro puede tener variaciones de 
hasta tres veces el valor nominal (K_. =3 
K.). 

e Por consiguiente, el parámetro de orden 
fraccional a del controlador FPI se utiliza 
para mejorar la robustez del sistema de 
control en lazo cerrado, en el sentido de 
la estabilidad, frente a variaciones de 
la ganancia K del Bocal; es decir, para 


garantizar que el margen de ganancia 


(mM) del sistema de control alcance su 
valor máximo. En la figura 6 se muestra 
la variación del margen de ganancia en 
función del parámetro a. De esta figura se 
observa que el valor máximo del margen 
de ganancia M, = 4.47 se alcanza con a. = 
0.77. Observar en esta figura que el margen 
de ganancia que se obtiene con a = 1 
(controlador Pl estándar) es M, =3.6, por lo 
que se logra una mejora de un 24%. 

e Frecuencia de cruce y tiempo de estable- 
cimiento: el tiempo de establecimiento 
de la planta nominal en lazo abierto es 
aproximadamente f%, = 3T,, + 1, = 3 000 
s. Se desea que el sistema de control en 
lazo cerrado sea dos veces más rápido, 
para ello se escoge la frecuencia de 
cruce 10, = 0.0033 rad/s, la cual de forma 
aproximada se corresponde con el tiempo 
de establecimiento t =5/1w, =1 500 s. Esta 
última expresión es una aproximación 
obtenida de la consideración general 
sobre la relación inversa que existe entre la 
frecuencia de cruce de ganancia y el tiempo 
de establecimiento (Ogata, 1993). 

e Margen de fase: se escoge un margen 
de fase q, = 60% el cual es un valor 
estándar para esta especificación. Para 
plantas de segundo orden, el coeficiente 
de amortiguamiento E se encuentra 

relacionado con el margen de fase 

mediante la expresión E = q, = 60/100 

si 0 < E < 0.6. Utilizando la expresión 


M, = exp [208 i=E ) la cual relaciona el 
máximo sobreimpulso M, y el coeficiente 


de amortiguamiento para plantas de se- 
gundo orden, se obtiene que el margen 
de fase seleccionado se corresponde con 
un M,= 10%, el cual es aceptable para 
esta clase de aplicación. Estas expresiones 
son utilizadas con frecuencia como 
aproximaciones razonables de plantas de 
alto orden (Ogata, 1993). 

Utilizando el par de especificaciones 


frecuenciales de diseño ((, , 1.) y el parámetro 
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Margen de ganancia 


Figura 6. Variación del margen de ganancia del sistema de control de la distribución de agua en el bocal, 


en función del parámetro al. 


de orden fraccional a, de las expresiones (8) y 
(9) se obtienen los valores de los parámetros 
de los controladores FPI y Pl, representados, 
respectivamente, mediante las siguientes ex- 
presiones: 


lA 
Ss 


Rep (s) = 0.026 (12) 


1+134.6 s 
S 


Rp, (5) =0.015 (13) 


Resultados de simulación 


Con el objeto de confirmar que el sistema de 
control de la distribución de agua en el Bocal 
presenta un mejor comportamiento, así como 
una mayor robustez con el controlador FPI 
que el controlador PI frente a variaciones 
de los parámetros dinámicos del Bocal o 
de perturbaciones externas, se realizaron 
diferentes simulaciones de dicho sistema 
con ambos controladores bajo las mismas 


condiciones de operación del canal. En la figura 
7 y en el cuadro 1 se muestran los resultados 
comparativos obtenidos de simulación, de 
donde se observa que cuando el Bocal presenta 
comportamiento nominal (régimen nominal de 
operación) ambos controladores se comportan 
de forma aproximadamente similar. Sin 
embargo, cuando se originan variaciones en 
el régimen de descarga del Bocal en el rango 
de operación [Q,... Q....] y como resultado la 
ganancia presenta un valor diferente al nominal 
(K € [0.01, 0.16]), la referencia de nivel de agua 
(320 cm) se alcanza mucho más rápido (menor 
tiempo de establecimiento) y con un menor 
pico máximo con el controlador FPI que con el 
controlador PI. Además, la atenuación al efecto 
nocivo de las perturbaciones es también mucho 
más rápida con el controlador FPI que con el Pl. 
Por consiguiente, el controlador FPI posibilita 
aumentar la operatividad de la distribución 
de agua en el Bocal y por ende disminuir las 
pérdidas por concepto de operación de los 
recursos hidráulicos disponibles. 
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322 Comportamiento 
nominal 

321 
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318 


Referencia / nivel de agua 


Perturbación 


Variaciones 


ed de la ganancia 


0 1 2 3 4 5 6 7 
Tiempo (s) x 10% 


Figura 7. Resultados comparativos simulados de la respuesta temporal del sistema de control de la distribución 
de agua en el bocal con controladores FPI y PL 


Cuadro 1. Resultados comparativos de simulación de la respuesta temporal del sistema de control de la distribución 
de agua en el bocal con controladores FPI y PL. 


as AK=-0.01 AK=0.01 As) =-0.3/8 
Controlador M, (%) t.(s) M, (%) t.(s) M, (%) t.(s) M, (%) t.(s) 
EPI 0.03 3 480 0.4 6 488 0.3 5120 0.3 2500 
PI 0.03 4 840 0.6 9320 0.7 6 844 0.3 4 000 


De esta forma se demuestra que con 
el controlador FPI se obtienen ventajas 


plementación práctica en tiempo real en un 
PLC Siemens SIMATIC S7-300, instalado en 


gua, vol. V, núm. 1, enero-febrero de 2014 
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significativas sobre el controlador PI cuando 
se requiere controlar la distribución de agua 
en un tramo de un canal principal de riego con 
parámetros dinámicos variables en el tiempo, 
comportándose de forma más robusta en todo 
el rango de variación de dichos parámetros. 


Resultados de implementación práctica 
del controlador FPI en el Bocal 


La viabilidad y robustez del controlador 
FPI diseñado se demostró mediante su im- 


la sala de control del Bocal, el cual forma 
parte del sistema de supervisión y control 
centralizado de la Confederación Hidrográ- 
fica del Ebro (CHE). Este PLC se caracteriza 
por ser compacto (modular) y muy potente, 
disponer de una alta velocidad de respuesta 
en tiempo real, presentar una alta conectivi- 
dad y ofrecer diferentes facilidades en el 
manejo del software y del hardware. Además, 
presenta alta inmunidad a las interferencias 
producidas por motores y puede trabajar 
en condiciones ambientales adversas, tales 
como altas temperaturas y elevada humedad. 
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El controlador PI (13) utilizado en el control 
del grupo de compuertas aguas arriba del Bocal 
(denominado regulador piloto) se encuentra 
implementado en el PLC Siemens SIMATIC S7- 
300 referenciado. Este controlador entrega su 
señal de control u,/(t) a otro controlador interno 
(denominado regulador de órgano), encargado 
de discretizar la señal de control u,/(kT) y 
obtener el posicionamiento (apertura real) del 
grupo de compuertas aguas arriba (situadas en 
la Casa de Compuertas). De esta forma se logra 
que las compuertas no estén continuamente en 
movimiento y se reduce el desgaste mecánico 
que un movimiento continuado trae aparejado. 
En la figura 8 se muestra el diagrama de 
bloques funcional del sistema de control de la 
distribución de agua en el Bocal, en el cual es 
posible observar la interacción existente entre 
el regulador piloto y el regulador de órgano. 

El controlador FPI diseñado (12) se 
implementó de forma paralela con el controla- 
dor PI (13) en el PLC SIMATIC S7-300 mediante 
un programa en AWL (lista de instrucciones). 
El operador de orden fraccional del controlador 
FPI (12) se implementó a través de un filtro 
discreto tipo FIR (Respuesta Impulso Finita) 
obtenido mediante su aproximación numérica 
utilizando la definición de  Grundwald- 
Letnikov (Podlubny, 1999a; Vinagre et al., 2000) 
y combinado con el principio de aproximación 
de memoria corta con serie de 500 términos (N 


= 500). El periodo de muestreo utilizado fue T 
= 60 s, por lo que el PLC dispone de suficiente 
tiempo para generar las señales de control. 

Los criterios que se siguieron para la im- 
plementación del controlador FPI (12) en el 
PLC SIMATIC S7-300 fueron los siguientes: a) 
estructura modular de las rutinas del código 
ejecutable, de forma tal que su funcionamiento 
disponga de la mayor autonomía posible; b) 
elevada fiabilidad, introduciendo en el PLC 
el código ejecutable estrictamente necesario 
para la implementación del controlador 
FPI de forma paralela con el controlador PI; 
c) simple y rápida conmutación entre los 
controladores FPI y PI en caso de detectarse 
algún tipo de fallo o de mal funcionamiento; 
y d) disponibilidad de las señales de control 
de ambos controladores, así como de la 
variable controlada (nivel de agua aguas abajo 
alejado) para la supervisión y comparación 
del funcionamiento de ambos controladores. 
El código ejecutable del controlador FPI (12) 
presenta las siguientes rutinas: 


a) Rutina de inicialización de datos internos. 
Se ejecuta sólo durante el arranque del 
sistema de control. Se ocupa de generar 
los vectores de pesos de la aproximación 
discreta de en 

función de los parámetros del controlador 


FPI (12) y del periodo de muestreo. 


Grundwald-Letnikov, 


Apertura inicial 
compuertas 


Consigna 
nivel canal + Regulador 


piloto (PI) 


Regulador de 
órgano 


Compuertas + 
canal + sensor 


Figura 8. Diagrama de bloques funcional del sistema de control de la distribución de agua en el bocal. 
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b) Rutina de selección del instante de 
muestreo. Se ejecuta de manera cíclica. Su 
función consiste en escoger uno de los dos 
controladores (FPI o PI) y en qué instante 
de tiempo se debe conectar. 

c) Rutina de cálculo de la señal de control. Se 
ejecuta a petición de la rutina anterior. Es la 
encargada de calcular la señal de control. 


Resultados experimentales 


Se muestran dos casos de resultados experi- 
mentales comparativos de funcionamiento 
en tiempo real del sistema de control de la 
distribución de agua con controladores FPI y 
PI en el Bocal. En el primer caso se consideró 
una disminución del valor de la señal de 
referencia de nivel de agua aguas abajo alejado 
(disminución de la descarga a través de las 
compuertas aguas arriba) y, en el segundo, 
un aumento de dicho valor (incremento de 
la descarga a través de las compuertas aguas 
arriba). En cada caso, se desarrolló primero un 
experimento con el controlador PI y luego se 
repitió el mismo experimento con el controlador 
FPI. Es necesario tener en cuenta que debido 
a que los experimentos con controladores FPI 
y PI se desarrollaron de forma separada en 
condiciones de operación en tiempo real es 
muy difícil mantener las mismas condiciones 
aguas arriba (variación de nivel en el río Ebro) 
durante ambos experimentos. De hecho, estas 
variaciones constituyen una de las principales 
perturbaciones del sistema de control en este 
tramo. Además, las variaciones de nivel en el 
río Ebro presentan carácter aleatorio (Pedregal 
et al., 2009) y no pueden ser controladas 
mediante nuestro sistema de control. 

En el primer caso se varió la señal de 
referencia de nivel de agua aguas abajo alejado 
del sistema de control de la distribución de 
agua en el Bocal desde 325 hasta 317 cm, lo 
cual representa una disminución de 8 cm, 
manteniéndose la compuerta aguas abajo en 
la misma posición. En este caso, la descarga 
a través de las compuertas aguas arriba 
disminuyó (Q,(t) <Q...) y, por consiguiente, el 
nivel de agua en el tramo también disminuyó. 


El gráfico superior de la figura 9 muestra los 
resultados comparativos dela variación denivel 
en el río Ebro (condiciones de perturbación), 
mientras que el gráfico inferior y el cuadro 2 
exhiben los resultados comparativos de la 
variación de nivel de agua aguas abajo alejado 
(variable controlada) en el Bocal durante los 
experimentos en tiempo real con controladores 
FPI y PI. Del gráfico superior de esta figura se 
observa que el nivel de agua en el río Ebro se 
mantuvo con una pequeña variación (Ay,(t) = 
+5 cm) durante cada uno de los experimentos 
con controladores FPI y PI, por lo que ambos 
experimentos se desarrollaron en condiciones 
similares de perturbación. Del gráfico inferior 
de esta figura y del cuadro 2 es posible observar 
que el nivel de agua en el Bocal durante ambos 
experimentos presentó una región de error 
de 2 cm en relación con el nivel de referencia. 
Sin embargo, el nivel de referencia se alcanza 
más rápido con el controlador FPI (52.5 min) 
que con el controlador PI (61.5 min), es decir, 
aproximadamente 9 min más rápido; a los 120 
min, el nivel de agua con el controlador FPI 
presenta cero error de estado estacionario; 
mientras que con el controlador Pl, esta 
variable todavía presenta un error de estado 
estacionario de 1 cm. 

En el segundo caso se varió la señal de 
referencia de nivel de agua aguas abajo alejado 
del sistema de control de la distribución de 
agua en el Bocal desde 310 cm hasta 320 cm, 
lo cual representa un aumento de 10 cm, 
manteniéndose la compuerta aguas abajo en 
la misma posición. En este caso, la descarga a 
través de las compuertas aguas arriba aumentó 
(Q,(t) > Q,...) y por consiguiente el nivel de 
agua en el Bocal también se incrementó. Los 
resultados experimentales comparativos de 
control de la distribución de agua en el Bocal 
con controladores PI y FPI se muestran en la 
figura 10 y en el cuadro 2. Del gráfico superior 
de esta figura (condiciones de perturbación) se 
observa un brusco descenso de nivel en el río 
Ebro durante el experimento con controlador 
FPI. Este descenso se produjo debido a la 
entrada en funcionamiento de una central 
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Figura 9. Resultados experimentales comparativos de funcionamiento en tiempo real del sistema de control de la distribución 


de agua en el bocal con controladores FPI y PI cuando la descarga disminuye. 


Cuadro 2. Resultados experimentales comparativos de funcionamiento en tiempo real del sistema de control de la distribución 
de agua con controladores FPI y PI en el bocal. 


Primer experimento Segundo experimento 
QUO < Quo OE 
Controlador t, (min) e. (cm) Ay, (cm) t.(s) e (cm) Ay, (cm) 
¡BN AS 0 ELO) 32 0 - 23 
PI 61.5 1 ED 48 0.8 E 


hidroeléctrica situada a la misma altura de 
la Casa de Compuertas del Bocal, pero en 
la margen contraria del río Ebro. Esta fuerte 
perturbación originó un descenso de nivel 
en el río Ebro desde los 385 hasta los 362 cm, 
es decir Ay,(t) = -23 cm (aproximadamente) 
en un tiempo inferior a los 2.5 min; mientras 
que durante el experimento con el controlador 
PI, el nivel en el río Ebro se mantuvo con 


una pequeña variación (Ay,(t) = +5 cm). Por 
consiguiente, el experimento con controlador 
FPI se desarrolló en condiciones de una 
mayor variación de la descarga a través de las 
compuertas aguas arriba del Bocal, así como 
de una mayor variación de los parámetros 
dinámicos de dicho tramo. 

En el gráfico inferior de la figura 10 y en 
el cuadro 2 se observa que el nivel de agua 
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Figura 10. Resultados experimentales comparativos de funcionamiento en tiempo real del sistema de control 


de la distribución de agua en el bocal con controladores FPI y PI cuando la descarga aumenta. 


aguas abajo alejado en el Bocal presentó una 
región de error de +1.0 cm en relación con el 
nivel de referencia durante los experimentos 
con ambos controladores, pero nuevamente 
el nivel de referencia se alcanzó más rápido 
(32 min) con el controlador FPI que con el 
controlador PI (48 min), es decir 16 min 
aproximadamente más rápido. Además, a los 
80.0 min, el nivel de agua en el Bocal con el 
controlador FPI presenta cero error de estado 
estacionario; mientras que con el controlador 
PI, esta variable todavía mantiene un error de 
estado estacionario de 0.8 cm. 

Estos resultados experimentales muestran 
que el controlador FPI exhibe un comporta- 
miento dinámico superior al controlador PI 
en el sentido de propiciar un menor tiempo 
de establecimiento de la respuesta temporal 


del sistema de control (mayor velocidad de 
respuesta) en aproximadamente un 25% y 
con una mayor exactitud en el control. Ello 
posibilita aumentar la operatividad de la 
distribución de agua en el canal, satisfacer en 
tiempo las demandas de los diferentes usuarios 
y reducir las pérdidas de agua por concepto 
de explotación en un 25% (Kovalenko, 1983; 
Litrico y Fromion, 2009). 

Por otro lado, se observa una mayor 
robustez del controlador FPI en relación con 
el PI debido a que la respuesta temporal del 
sistema de control con controlador FPI presenta 
un mejor comportamiento que la obtenida con 
el controlador PI bajo condiciones de mayores 
variaciones de los parámetros dinámicos del 
Bocal, así como de perturbaciones externas 
más agresivas. 
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De esta forma se demuestra con resul- 
tados experimentales comparativos de funcio- 
namiento en tiempo real del sistema de control 
de la distribución de agua en el Bocal con 
controladores FPI y PI la conveniencia de la 
implementación de controladores de orden 
fraccional cuando se requiere el control efectivo 
de canales principales de riego con parámetros 
dinámicos variables en el tiempo. 


Conclusiones 


Se presentó el diseño e implementación prác- 
tica en tiempo real de un controlador FPI en 
un tramo de un canal principal de riego en 
explotación, el cual se comporta de forma más 
robusta que otro controlador PI equivalente 
(en el sentido de que ambos exhiben el mismo 
comportamiento dinámico en lazo cerrado 
para las especificaciones de la planta nominal) 
frente a grandes variaciones de los parámetros 
dinámicos, así como ruidos de alta frecuencia 
de la planta. Este resultado presenta una 
especial relevancia en canales principales de 
riego, cuyos parámetros dinámicos exhiben 
una amplia variación cuando se originan 
cambios en los regímenes de descarga en el 
rango de operación [Q,.., OQ... ]. 

Se desarrolló la implementación práctica 
del controlador FPI diseñado en un PLC 
SIMATIC 57-300, con base en criterios de 
simplicidad, elevada autonomía y fiabilidad 
del código ejecutable del controlador. 

Se desarrollaron experimentos evaluativos 
del funcionamiento en tiempo real del sistema 
de control, con controladores FPI y PI de la 
distribución de agua en un tramo de un canal 
principal de riego en explotación, caracterizado 
por presentar amplias variaciones en sus 
regímenes de descarga. Estos experimentos 
mostraron que cuando se originan variaciones 
en el régimen de descarga del canal (Q,(t) 
+ Q.,), con el controlador FPI se obtiene 
un menor tiempo de establecimiento de 
la respuesta temporal (mayor velocidad 
de respuesta) del sistema de control en 
aproximadamente un 25% y con una mayor 


exactitud en el control, lo cual implica una 
mayor operatividad de la distribución de agua 
en el canal, satisfacer de forma más rápida las 
demandas de los diferentes usuarios y reducir 
de forma significativa las pérdidas de los 
recursos hidráulicos disponibles por concepto 
de operación en aproximadamente 25%. 

Resulta importante destacar que en este 
trabajo se presentan las primeras experiencias 
relacionadas con la implementación práctica 
y evaluación comparativa de un controlador 
robusto de orden fraccional FP] en un tramo de 
un canal principal de riego en explotación. El 
siguiente objetivo de nuestras investigaciones 
consiste en generalizar los resultados obtenidos 
en este trabajo en el control integral de un 
canal principal de riego conformado por varios 
tramos con fuerte interacción entre los tramos. 

Los beneficios que se obtienen con el 
aumento de la robustez y efectividad en 
el control de la distribución de agua en un 
canal principal de riego se revierten de forma 
inmediata en una mejor administración y 
uso de los recursos hidráulicos disponibles, 
así como en una mayor protección del medio 
ambiente, al disminuir las pérdidas de agua 
por concepto de operación. 
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